Appearance of room temperature ferromagnetism in Cu-doped Ti02-5 films. 
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We report here the unexpected observation of significant room temperature ferromagnetism in a 
semiconductor doped with nonmagnetic impurities, Cu-doped Ti02 thin films grown by Pulsed Laser 
Deposition. The magnetic moment, calculated from the magnetization curves, resulted surprisingly 
large, about 1.5 fiB per Cu atom. A large magnetic moment was also obtained from ab initio 
calculations , but only if an oxygen vacancy in the nearest- neighbor shell of Cu was present. This 
result suggests that the role of oxygen vacancies is crucial for the appearance of ferromagnetism. 
The calculations also predict that Cu doping favors the formation of oxygen vacancies. 



Integrating spin functionality into otherwise nonmag- 
netic materials has become a highly desirable goal in 
the last years. For example, dilute magnetic impurities 
in semiconductors produce novel materials appealing for 
spintronics (see, e.g., refs. [l|,0, Hand references therein). 
This is a rapidly developing research area, in which the 
electron spin plays an important role in addition to the 
usual charge degree of freedom. For their practical ap- 
plications, ferromagnetic semiconductors are required to 
have a high Curie temperature (Tc). While most of the 
dilute magnetic semiconductors (DMS) like Mn-doped 
GaAs Q have a Tc much lower than room temperature, 
room-temperature ferromagnetism has been observed in 
some Mn-doped compounds such as ZnO [f|. Recently, 
Co-doped Ti02 thin films with the anatase structure were 
reported to be ferromagnetic even above 400 K with a 
magnetic moment of 0.32 \ib per Co atom jf|. These 
results have motivated intensive experimental 0, S, 0| 
and theoretical 13, U, 13, 13, 14 1 studies on the struc- 
tural and electronic properties of these materials. How- 
ever, many questions remain open regarding the underly- 
ing microscopic mechanism of long-range magnetic order. 
Carrier-induced interaction between the magnetic atoms 
was first suggested as the important in gred ient underly- 
ing ferromagnetism in III-V based DMS [15| . Subsequent 
reports have raised concerns about the initially suggested 
intrinsic nature of ferromagnetism in these materials, due 
to the possibility of ferromagnetic metal clustering under 
different growth conditions [EL f]~6j] . Furthermore, it has 
been suggested that the strong interaction between tran- 



sition metal impurities and oxygen vacancies is crucial for 
the existence of high Curie temperature ferromagnetism 
18, [ljjj]. We believe our results shed new light on 
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this problem, suggesting a specific role for the oxygen 
vacancies and against the clustering hypothesis. 

The original purpose of this work was to find out the 
role of oxygen deficiency in the origin and significance 
of ferromagnetism in doped Ti02 thin films. Thin films 
are used because they have a much greater surface or in- 
terface to volume ratio than single-crystals or polycrys- 
talline ceramics, so the role of defects may be enhanced. 
On the other hand, we had previously found in Fe-dopcd 
TiC>2 thin films analyzed by X-Ray Absorption Near- 
Edge and Extended X-Ray Absorption Fine Structure 
Spectroscopy that the magnetic signal increased when 
the number of oxygen vacancies around the impurity sites 
increased [2(|. In order to produce a similar concentra- 
tion of vacancies in Ti02 thin films as in Fe-doped sam- 
ples but without the magnetic ions, we studied Cu-doped 
Ti02 films assuming that no magnetic signal would come 
from the dopant. Surprisingly, significant room tempera- 
ture magnetic behavior, so strong to give a magnetization 
equivalent to 1.5 ^b/Cu, was found. 

Thin films of approximately 10 at.% Cu-doped Ti02_5 
were deposited on LaA10 3 (001) substrate (LAO) by 
Pulsed Laser Deposition (PLD) using a Nd:YAG laser 
operating at 266 nm. The Ti02:Cu target was prepared 
from high purity Ti02 and metallic Cu powders in stoi- 
chiometric quantities. The powders were mixed for three 
minutes using a ball-mill, then uniaxially pressed (200 
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FIG. 1: X-Ray Diffraction pattern for Cu-doped T1O2-S film 
deposited by PLD on a LaAlOs substrate. 



MPa) into a disk, and finally sintered. The substrate 
temperature, laser energy density, oxygen pressure, and 
pulse repetition rate were 800 °C, 2 J/cm 2 , 20 Pa, and 
10 Hz, respectively. After deposition the film was cooled 
down to room-temperature in 2 hs under 20 Pa oxy- 
gen atmosphere. The composition was determined by 
Energy Dispersive X-Ray Analysis (EDX), and no con- 
taminants were found within the accuracy of the method 
(< I wt. %). The crystalline structure was studied by 
X-Ray Diffraction (XRD). Our films were transparent 
and strongly textured, showing only the (001) reflections 
of the anatase structure, although rutile reflections were 
also present but with much less intensity (Fig. [T|). This 
strong difference between intensities does not directly re- 
flect the relative amount of both phases, as our samples 
consist of a highly oriented anatase film [only (001) re- 
flections with very large intensities] and polycrystallinc 
rutile embedded in it [(110) and (101) reflections of low 
intensities]. The relative amount of the anatase and ru- 
tile polymorphs depends on the choice of substrate and 
on the growth conditions. The lattice mismatch between 
anatase Ti0 2 (001) and LaA10 3 (001) is only 0.26%, 
much less than for rutile, so anatase is favored when de- 
position is performed on (001) LAO substrate. However, 
it has been shown [2l| that a slow growth rate is required 
to produce anatase single crystal films of high crystallo- 
graphic quality, while growing at a higher rate results in 
polycrystallinc films containing both anatase and rutile 
phases. Chambers et al. 2l| also found that the rutile 
nanocrystals nucleate within the evolving anatase film in 
orientations that maximize the lattice match between the 
two phases. 

The measurements of magnetization M as a function 
of the applied magnetic field H were performed with a 
commercial Vibrating Sample Magnetometer Lake Shore 
7407 at room temperature, with the external field ap- 
plied parallel and perpendicular to the plane of the film. 



FIG. 2: Room temperature hysteresis curve of the Cu-doped 
T1O2-S film (points) in the sample as cast and after a thermal 
treatment (TT) at 800 °C for 30 min in a oxygen-rich atmo- 
sphere. Solid lines are the best least-squares fits of Eq. [1] to 
the data. The inset compares the hysteresis loops measured 
with the external magnetic field applied parallel (PA) and 
perpendicular (PP) to the film plane for the as cast sample. 



After subtraction of the diamagnctic contribution due to 
LaAlOs we obtained the result depicted in Fig. [2j where 
significant room temperature magnetization is displayed. 
In order to quantify the magnetic parameters (saturation 
magnetization M s , intrinsic coercivity H c , and remanent 
magnetization M r ) we propose the following fitting func- 
tion for the demagnetization data: 



M = Af s [(2/7r)[arctan((i? + H c )/H c ) tan(7rS/2)]] + \H, 

(1) 

where S = M r /M s . The first term is the usual func- 
tion used to represent a ferromagnetic hysteresis curve 
[23I ] and the second is a linear component representing 
a possible paramagnetic contribution. The parameters 
obtained are: 44 8 emu/cm 3 , 2.54 emu/cm 3 , and 353 Oe 
for M s , M r , and H c , respectively. From M s , we estimate 
an atomic moment of 1.5 /kb/Cu, assuming 10 at. % of 
Cu and a 1000 A film thickness. 

For the interpretation of this unexpected experimental 
result we performed ah initio calculations on doped Ti02, 
with and without oxygen vacancies, using the supercell 
(SC) method. As a first approximation we employed a 
small SC (T14O8) consisting of two unit cells of the ru- 
tile structure stacked along the c-axis. One Ti atom was 
replaced by Cu, giving a rather large impurity concen- 
tration compared with the experimental one. However, 
the appearance of a magnetic moment in some cases was 
also quite surprising and interesting to relate with the 
experiments. The rutile structure was chosen for the cal- 
culations because it gives a more uniform distribution of 
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defects than the anatase structure when the SC is small. 
This leads to a more realistic approximation of the struc- 
tural distortions induced by the impurity and the vacancy 
in the host lattice. The calculations were performed 
with the ab initio full-potential linearized-augmented- 
plane-wave method (FP-LAPW) in a scalar relativistic 
version, as implemented in the Wien2K code [24|. We 
used the local spin density approximation (LSDA) and 
the exchange-correlation potential given by Perdew and 
Wang (ill. The cutoff criterion was R mt KMAX=7, be- 
ing R mt the smallest muffin tin radius and Kmax the 
largest wave number of the basis set. The number of k- 
points was increased until convergence was reached. We 
also employed, specially for structural minimizations, the 
SIESTA code [26| that uses a linear combination of nu- 
merical real-space atomic orbitals as basis set and norm- 
conserving pseudopotentials. We proved using Wien2K 
that the relaxed structures predicted by SIESTA have 
lower energies than the unrelaxed ones and that both 
methods give similar magnetic moments. 

The SC used for these calculations is tetragonal 
(a=6=4.5845i A, c'=2c=5.9066 2 AQ) and contains 4 
Ti atoms at [(0,0,0); (1/2,1/2,1/4); R] and 8 O atoms at 
Uu, u, 0); (1/2 + u, 1/2 - u, 1/4); R] with u=0.3049i 
[22|. The notation R indicates that the group of coor- 
dinates should be repeated adding (0, 0, 1/2). In this 
structure, the Ti atoms are surrounded by a slightly 
distorted oxygen octahedron with a rectangular basal 
plane of four oxygen atoms with distances to Ti (1.94 
A) slightly shorter than those at the vertex (1.98 A). 

For the case of Ti 3 CuOs, the system without oxygen 
vacancies, we found that the presence of Cu induces a 
small distortion of the host lattice, mainly in the six oxy- 
gen nearest-neighbors (Oat at) of the Cu impurity. The 
Cu-OjVAr distances are enlarged from 1.94 A and 1.98 A 
to 1.96 A and 2.01 A, respectively. However, no mag- 
netism was found, which suggests that the experimen- 
tally observed magnetic moment cannot be due only to 
the presence of Cu impurities. 

In order to analyze the effect of oxygen vacancies in 
doped Ti02 we first studied the system T14O7, taking 
into account the structural distortions produced by the 
vacancies. These are much larger than those found in 
Ti 3 CuOs, as some of the Ti atoms are displaced as much 
as 0.2 A. We found that for this system there are two 
possible degenerate solutions, one is non magnetic and 
the other has a magnetic moment of 1.1 /is per SC, due 
to polarization of the Ti atoms. It is worth mentioning 
that the magnetic phase appears only when structural 
relaxation is taken into account. This result anticipates 
that oxygen vacancies could play an important role in the 
magnetic behavior observed in the Ti02:Cu films. The 
question that arises is why the magnetic phase is not ob- 
served in bulk rutile TiC>2 without impurities. A possible 
explanation for this is that the concentration of intrinsic 
vacancies formed in equilibrium is not enough to generate 



magnetism at the temperature considered. Calculations 
with a larger unit cell and only one vacancy, both in ru- 
tile and anatase structures, were studied and a smaller 
magnetic moment was obtained for the rutile structure 
and none for the anatase structure. So it seems that the 
appearance of a magnetic moment may depend on the va- 
cancy concentration and on the structure. An experimen- 
tal observation of magnetism in an undoped nonmagnetic 
oxide has been recently observed experimentally and has 
been attributed to vacancies also 27 1. 

Finally, we studied the oxygen deficient Cu-doped 
TiC-2 systems (Ti 3 Cu07). We considered the three pos- 
sibilities for removing an oxygen atom from the SC: 
(i) from the Ti-contained octahedron, (ii) from the Cu- 
contained octahedron removing a basal O atom, and (iii) 
from the Cu-contained octahedron removing a vertex O 
atom. Comparing the total energies of the 3 cases we 
found that the oxygen vacancies near Cu are more stable 
than those near the Ti site. A similar result was found 
by Weng et al. [l3| for Co in oxygen-deficient Ti02- 

The calculations also predicts that the energy required 
to form an oxygen vacancy decreases from 10 eV in the 
undoped system to 4 eV in the case it is near a Cu impu- 
rity. Also, the energy required to replace a Ti atom by 
a Cu atom is reduced by 5 eV when there is an oxygen 
vacancy in the SC. We therefore predict that Cu doped 
systems will have more vacancies than the undoped ones, 
and that oxygen vacancies will tend to be close to the im- 
purities. We also found that the energy required to form 
two vacancies near a Cu impurity is larger than that nec- 
essary to create one vacancy in the first shell of neighbors 
of two different Cu impurities. For this reason, we discard 
the formation of vacancy clusters around Cu atoms. 

Concerning the magnetic properties of the system 
Ti 3 Cu07, we found that the magnetic solution has lower 
energy than the nonmagnetic one, with a magnetic mo- 
ment of 1.0 fiB per SC, due to polarization of both Cu 
and Oat at atoms. The value of the magnetic moment is 
nearly independent of vacancy location and of structural 
relaxation. This result indicates that there is magnetism 
in Cu-doped Ti02, but as we have considered only one 
SC size (and therefore only one impurity concentration) 
we cannot compare its value with the experimental mag- 
netic moment per Cu atom. To do this we should average 
over three different concentrations: that of Cu, that of 
vacancies, and the relative amount of rutile and anatase 
phases. As a first step in this direction we performed a 
few calculations with a larger SC containing 48 atoms. 
The dimensions of this SC are a'=6'=2a=26=9.169 A, 
c'=2c=5.907 A. In particular, for the case of one Cu atom 
with a neighbor oxygen vacancy in the SC we obtained 
a magnetic moment of approximately 1.0 /is, which val- 
idates the results obtained for the smaller SC. A similar 
result was found for the anatase phase. 

Fig. [3] shows a calculated density of states of the sys- 
tem Ti 3 Cu07, which is small gap semiconductor or a 
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a magnetic moment of 1.0 \ib for a supcrcell containing 
one Cu impurity and a neighbor oxygen vacancy, but no 
magnetic moment if the oxygen vacancy is absent. The 
calculations also predict a lower vacancy formation en- 
ergy when there are Cu impurities. Therefore, it seems 
that oxygen vacancies play a crucial role for the appear- 
ance of magnetism and that one effect of doping is to 
increase their number. Both experimental and theoreti- 
cal work is now in progress. 

This work was partially supported by ANPCyT, 
PICT98 03-03727, CONICET, Fund. Antorchas, Ar- 
gentina, and TWAS, Italy, RGA 97-057. We thank R. 
Weht for his assistance with the SIESTA code. 



FIG. 3: Density of states of TigCuOy, with the oxygen va- 
cancy near Cu. Positive values are for spin up and the neg- 
ative ones for spin down. Energies are referred to the Fermi 
level, Ef. The minority spin feature inside the gap but close 
to Ef is 50% of Cu 3d and 50% of oxygen and Ti character. 



semimetal. However, if the number of vacancies were 
slightly larger than that of the Cu impurities, the Fermi 
level would move into the minority spin feature close to 
it. In this case the system would be magnetic and at the 
same time half-metallic, with carriers of only one spin 
type; this could be very interesting for spintronics. A 
similar idea, with excess vacancies responsible for ferro- 
magnetism was proposed recently [28j |. 

In the previous calculations a ferromagnetic alignment 
of the Cu spins was assumed, but we also verified that 
the pairs Cu-Oat^ vacancy do not couple antifcrromag- 
netically. This last calculation was performed for a few 
particular cases, as the number of possible distributions 
of impurities and vacancies is large and so is the supercell 
required for that calculation. 

In order to verify the crucial role played by the oxygen 
vacancies in the origin of ferromagnetism, we performed 
a "reverse experiment": in order to reduce the oxygen 
vacancies, the Ti02 film was heated at 800 °C for 30 min 
in a oxygen-rich atmosphere. While the XRD spectra 
show no changes, the measurements of M as a function 
of H after the thermal treatment show a drastic reduction 
of the magnetic signal (see Fig. [2]), in agreement with the 
ab initio predictions. 

In conclusion, we have observed an unexpected signif- 
icant room temperature magnetic behavior in Cu-doped 
Ti02-<5 films, so strong to give 1.5 [Lb per Cu atom. This 
experimental result indicates that magnetic ions are not 
essential to obtain this effect and also that is not due to 
impurity clustering. Ab initio calculations on bulk Cu- 
doped Ti02, with and without oxygen vacancies, predict 
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